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Abstract 


The aim of this article is to present a discrete model for histogram shaping. This is 
an important image transformation with several practical applications. The model that is 
proposed is based on a generalization of the inferior part function. Finally, an algorithm 
based on this model is developed. 
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1 Introduction 


Histogram equalization or histogram flattening is one of the most important nonlinear point 
operations. This transformation aims to distribute uniformly the gray levels of the input image 
such that the histogram of the output image is flat. Histogram equalization has been studied 
for many years (see [1], [3], [4]) and many practical applications have been proposed so far. A 
direct generalization of this transformation is represented by histogram shaping or histogram 
specification (see [1], [3]). The idea of histogram shaping is to transform the input image 
into one which has histogram of a specific shape. Obviously, when the output shape is flat, 
histogram equalization is obtained. Both histogram equalization and histogram shaping have 
become classical image transformations, therefore it has been quite difficult to find the initial 
reference source. One of the earliest references about is (2]. 

The mathematical model of histogram equalization and shaping is based on stochastic ap- 
proach. Let us consider that the input digital image is f = (f; g: t=1,2,...,.n; j =1,2,...,.m) 
where 

l<fig<G 
represents the gray value of pixel (i, ). The probability or frequency of gray level k € 1,...,G is 
defined by 


p(k) = 2 ==) ey ee = FeawG, (1) 


where #{(¢, 9): f(i,7) =k} gives the number of pixel with the gray level equal to k. Based on 
these probabilities, the digital image f can be considered a discrete random variable 


mo 2 a1 
= (Joan py(2) - py(G) } @) 


for which ey p;(k) = 1. Recall that the cumulative probability distribution of pf is 
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k 
Pr: {1,2,..., G} — [0,1], Pr(k) = S~ p,(i). (3) 
1=0 
A more productive approach is to consider the digital image as a continuous random variable py : 
[0, 00) — [0,1] such that {5° ps(x)da = 1. In this case the cumulative probability distribution is 


P; : [0,00) — [0,1), P;(z) = i, ” py (t)dt. 


Based on this continuous model the histogram shaping transformation can be defined more 
easily. Consider that the input digital image f is transformed such that the histogram of the 
output image g has a shape given by the cumulative distribution Q : [(0,c0) — (0,1). The 
equation that gives histogram shaping is [3] 


9=Q-*(P(f)). (4) 


The main inconvenience arising form Equation(4) is represented by the inverse function Q-?, 
Firstly, because the calculation of Q—! might not be easier even for simple shapes. Secondly, we 

_cannot define Q! for the discrete case therefore it would be difficult to apply (4) to a discrete 
computation. Perhaps, this is the real reason for seeing no discrete models for histogram shaping. 
In the following we will propose a discrete model for this transformation. 


2 The Superior Smarandache f-Part 


In order to propose an equation for the discrete case, we have to find a substitute for Q-!. This 
is given by the Superior Smarandache f-Part , which represents a direct generalization of the 
classical ceiling function. Smarandache proposed [5] a generalization of the ceiling function as 
following. Consider that f : Z — R an increasing function such that limn—.—oof (n) = —oo and 
litn—oof(n) = oc. The Superior Smarandache f-Part. associated with f is fll: R + Z defined 
by 
fic) =k @ f(k-1)<2< f(h). 

Smarandache studied this function in relation to some functions of N umber Theory and proposed 
several conjectures on them [6]. Tabirca also studied the Superior Smarandache f-Part [7] when 
f(n) = So # and proposed equations for f" when a = 0,1,2. Tabirca also applied this 
function to static parallel loop scheduling [8]. 

Now, we propose a version of the Superior Smarandache f-Part for our discrete case. Consider 
that f : {1,2,...,G} — (0, 1] is an increasing function such that f(G) = 1. We also consider 
that this function is extended to 0 with f(0) =0. The Superior Smarandache f-Part associated 
with f is fl: 0,1] 4 {1,...,G} defined by 


fla) =k @ f(k—1)<2< f(k), Vz € [0,1]. (5) 


This function is also extended in 0 by f1(0) = 0. 
Some properties of the function f! are proposed in the following. 


Theorem 1 
f" (f(k)) =k, Vk € {1,2,...,G}. (6) 
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Proof The proof is based on the definition of f! and on the double inequality 
f(R-1) < f(k) < F(R). 


r 


Theorem 2 
a <f(fl(@)) <o+sup(F(k +1) — F(R), Ve € (0,11. (7) 


Proof 

Let us denote k = fl(x). The definition of fl provides f(k — 1) < 2 < f(k). From this 
equation, it directly follows that x < f(fll(2)). . 
The second part of Equation (7) comes from the following implication: 


i(k) < f()+2-fk-1) > 
§ (#1@)) < 2+ sup (F(k) — Fk). 


& 
Based on these properties, the histogram shaping model of the discrete case is proposed. 


3 Histogram Shaping for the Discrete Case 


Consider that the input image f = (f:;: i=1,2,...,n; j = 1,2,...,m) is transformed into the 
output image g = (9:3: i1=1,2,...,n; 7 =1,2,...,m) such that the histogram of g has a certain 
shape. Let us presume that the shape of the output histogram is given by the discrete random 
variable 


1 2 . G 
m= (0) Pa(2) - pr(G) } ®) 


where )-¢_,p,(k) = 1. 
The general equation of histogram shaping is similar with Equation (4) but pi is used in 
place of P,'. Let consider that the equation of image g is 


g(t, 5) = PE (Ps(F4,9))), V G5) € (1, -.0} x {1,2,...m}. (9) 


We prove that the cumulative probability distribution of g is very close to the cumulative 
probability distribution of h. 


Theorem 3 
P,(k) = Py (PH (Pa ())) , Wk € {1,2,...,G}. (10) 


Proof The proof is given by the following transformations: 


k 
Py(k) = Prig(i,j) <] = SoPr [Pl (P¢@9)) =]] = 
I=1. 
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k 
= 5° Pr[Pa(l—1) < Py(f(i, 9) < Pa(D] = 
l=1 


k 
= 0Pr [Py (Pall -1)) < PF (Py(FG,9))) < PP PAO)] = 
i=1 


=> (Pe[s6.0) < PE PO)] ~Pr [6.9 < PPG i) = 


I=1 


= >> (P; (P} (PD) - P; (P(A -1)))) = 


I=] 
= Py (PF (Pa())) — Py (PP (PA(0))) 
Since P; (P} (P,(0))) = 0 we find that P,(k) = Py (P} (Pa(k))) holds. 
From Theorems 2 and 3 the following equation is directly obtained. 


Pa(k) < Po(k) < Pa(k) + sup (P¢(j + 1) — Py(7)), VK € {0,1,...,@ — 1}. (11) 


Equation (11) provides an estimation of the gap between the quantities P;,(k) and P,(k). When 
sup; (P;(j +1) — Py(j)) is smaller these two quantities are very close. Although Equation (11) 
does not give a perfect equality we can say that the histogram of the image g has the shape very 
close to h. 

The algorithm based on this model firstly finds the functions Py, Pp, and Pi. Secondly, 
Equation (9) is applied to obtain the value g(i, j) for each pair (i,j). A full description of this 
algorithm is presented below. 


Inputs: 

n,m - the image sizes. 

f=(f[4,j]: i=1,...,n; j=i,...,m) - the input image. 

p_h=(p_h[i]: i=1,...,G) - the desired shape. 
Output: 

e=(gli,j]: i=1,...,n; j=1,...,m) - the input image. 
double P_h(int k){ 

double s=0; 

if (k<=0 [| k>G) return 0; 

for(int i=1;i<=k;i++)s=s+p_h[i]; 


return s; 

} 

int P_h_Inv(double x){ 
int k; 


“if (x<=0 || x>1) return 0; 
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Figure 1: Lenna’s picture 1. 


for(k=0;P_h(k)<=x;k++); 
return k-1; 


double P_f(int k){ 
double s=0, p_f[G]; 
if (k<=0 || k>G) return 0; 
for(int i=1;i<=G;i++)p £[i]=0; 
for(int i=1;i<=n;it++)for(int j=1;j<=m; j++)p_f[tli,j]]=p_t (tli, jl]+1./(nsm); 
for(int i=1;i<=k;i++)s=s+p_h[i]; 
return s; 
} 
int []0) Hist_Shaping(int n, int m, int fOO)L 
int g{n] fn]; 
for(int i=1;i<=n;i++)for(int j=1; j<=m; j++) 
gli,j] = P_h_Inv(P_f(f£[i,j])); 
return g; 


In order to show that the algorithm performs well we consider an example presented in [1]. 
Histogram shaping can be used to compare two images of the same scene, which have been taken 
under different lighting conditions. When the histogram of the first image is shaped to match 
in the histogram of the second image, the lighting effects might be eliminated. 

Consider that we have the images presented Figures 1 and 2. They are two different Lenna’s 
images where the second one has a poor lighting. Each image also contains the histogram for 
the red channel. The histogram shaping algorithm was applied to transform the second image 
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Figure 2: Lenna’s picture 2. 


according to the histogram of the first image. Figure 3 shows the resulting image which is the 
same as the first image. Moreover, the histograms of the first and third images are very alike 
with similar positions for peaks and valleys. 


4 Conclusions 


This article has introduced a discrete model for the histogram shaping transformation. The 
model that has been proposed uses the Smarandache ceiling function and is based on the equation 
g= Pil (P;(f)). A example has been also presented in order to prove that the method is viable. 
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Figure 3: Lenna’s Picture after Histogram Shaping. 
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